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ABSTRACT 

We study the effects of galaxy formation on the Sunyaev-Zel’dovich effect (SZE) observable-mass 
relations using high-resolution cosmological simulations. The simulations of eleven individual clusters 
spanning a decade in mass are performed with the shock-capturing Eulerian adaptive mesh refinement 
N-body+gasdynamics ART code. To assess the impact of galaxy formation, we compare two sets of 
simulations performed in an adiabatic regime (without galaxy formation) and those with several 
physical processes critical to various aspects of galaxy formation: radiative cooling, star formation, 
stellar feedback and metal enrichment. We show that a SZE signal integrated to a sufficiently large 
fraction of cluster volume correlates strongly with its enclosed mass, independent of details of gas 
physics and dynamical state of the cluster. The slope and redshift evolution of the SZE flux-mass 
relation are also insensitive to processes of galaxy formation and are well characterized by a simple 
self-similar cluster model. Its normalization, on the other hand, is significantly affected by gas cooling 
and associated star formation. Our simulations show that inclusion of these processes suppresses the 
normalization by « 30 — 40%. The effect is due to a decrease in gas mass fraction, which is offset 
slightly by an increase in gas mass-weighted temperature. Gas cooling and star formation also cause an 
increase in total mass and modify the normalization by a few percent. Finally, we compare the results 
of our simulations to recent observations of the SZE scaling relations obtained using 36 OVRO/BIMA 
SZE+Chandra X-ray observations. The comparison highlights the importance of galaxy formation in 
theoretical modelling of clusters and shows that current generation of simulations produce clusters 
with gross properties quite similar to their observed counterparts. 

Subject headings: cosmology: theory-clusters: formation- methods: numerical 


1. INTRODUCTION 

The Sunyaev-Zel’dovich effect (SZE) is a potentially 
powerful observational tool for cosmology. It is a small 
distortion in the cosmic microwave background (CMB) 
spectrum caused by scattering of CMB photons off a 
distribution of high energy e lectrons in dense struc tures 
such as clusters of galaxies (iSunvaev fc ZeldovichHloYfi 
dil). This effect has the unique property that its sig¬ 
nal is independent of redshift, makin g it partic ularly well 
suited for deep clu ster surveys (e.g./ lHolder et al ■ '2000: 
iWeller et al .11200211 . The next generation of SZE instru¬ 
ments, such as the South Pole Telescope (SPT) and the 
Atacama Cosmology Telescope (ACT), should be capa¬ 
ble of mapping a fairly large portion of the sky and 
finding a large number (> 10 4 ) of clusters out to high- 
redshift. Such large and homogeneous sample of galaxy 
clusters will enable direct and precise measurements of 
their number density as a function of redshift, and the 
expected survey yields will be sufficient to provide one 
of the most powerful constraints on the nature of da rk 
energy llWang k Steinha.rdtll 1 9981 lHaiman et al .11200 H) . 

To realize the full statistical power of the upcoming 
SZE surveys, however, systematic uncertainties would 
have to be controlled at a level comparable to statistical 
uncertainties. One of the main sources of the systematic 
uncertainties lie in the relation between the SZE observ¬ 
able and cluster mass as a function of redshift. Making 
this connection is important because the cluster mass is 


not directly observable. For the future SZE surveys, the 
requirement is to control systematic uncertainties in the 
SZE observab le-mass relations to be tter than ~5% at all 
redshift (e.g.. ICarlstrom et al .1120(14l . 

This poses a serious challenge to both observers 
and theorists. To date, observational studies of the 
SZE scaling relations have been performed using two 
largest datasets of SZE measurements obtained by the 
OVRO/BIMA cm-wave imaging experime nt (ICoora.vl 
119991: IMcCarthv et al.l 120031: lLaRoauel||2005l) and multi¬ 
band SuZIE experiment llBenson et al.l[2004[) along with 
X-ray observations. Analyzing a sample of 36 clusters 
observed with both the OVRO /BIMA S Z E ima ging and 
Chandra X-ray observations, ILa.Romifi ll2005lf showed 
that there are tight correlations between the observed 
SZE flux and X-ray temperature and cluster masses. The 
observed regularity of the SZE effect of clusters is encour¬ 
aging news, but further progress is clearly needed for the 
future SZE surveys. The observational situation is ex¬ 
pected to improve rapidly with the advent of a number 
of dedicated SZE survey instruments, which will dramati¬ 
cally increase the sample size and the number of low-mass 
clusters. 

On a theoretical side, a number of groups have stud¬ 
ied the SZE scaling relations using semi-ana lytic mod¬ 
els llVerde et al.l l2004 IMcCarthv et all I2003D and cos¬ 
mological simulation s iMgJzlgi^^ ; |White_C^l|j2^2 
da Silva et all 12004 iDiafarioetal 120051 iMotle^thl 
200511 . Motivated by results of X-ray observations in a 
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past decade (see IVnitll200'5i. for a review and references 
therein), recent studies have focused on studying the ef¬ 
fects of non-gravitational physical processes, including 
gas cooling, star formation and energy feedback, on the 
SZE scaling relations. One of the main results is that 
a SZE signal integrated to a sufficiently large fraction 
of cluster volume is an extremely good proxy for its en¬ 
closed mass, independent of details of gas physics an d 
dynamical state of a cluster (see e.g., IMotl et aHl2005 (). 
The slope and the redshift evolution of the SZE scaling 
relation also appe ar to be insensitive t o details of cluster 
gas physics (e.g.. Ida Silva et all 12004 1. While these re¬ 
sults are encouraging for cosmological applications, these 
previous studies have focused on simulating a large num¬ 
ber of clusters, and the resolution was inevitably limited 
to capture relevant cluster physics. As such, the impact 
of galaxy formation on normalization of the SZE scaling 
relations has not yet converged among different simula¬ 
tions. It is therefore important to push theoretical mod¬ 
elling of SZE scaling relations and check previous results 
using higher-resolution cluster simulations. 

In this paper, we present such study using high- 
resolution cosmological simulations of cluster formation. 
Although the statistic is limited, our cluster sample spans 
over a decade in mass and provide a good leverage on 
scaling relations. The mass resolution of our simulations 
is more than an order of magnitude higher than that in 
previous studies. This work is therefore complimentary 
to the previous studies in a literature. Using these sim¬ 
ulations, we study the impact of gas cooling and star 
formation on the SZE scaling relations, including their 
normalization, slope and redshift evolution. To test the 
results of our simulations, we compare our results to re¬ 
cent observations of the SZE scaling relations based on a 
sample of 36 clusters obtained using the OVR O+BIMA 
SZE and Chandra X-ray telescopes llLaRoauelll2005 f. 

The paper is organized as follows. In § [2] we define 
observational quantities and present relevant scaling laws 
predicted by a self-similar cluster model. We describe 
simulations presented in this paper in § 0 and present 
results and comparisons to previous studies and recent 
observational results in § Q] Finally, in § 0 we discuss 
our conclusions and their implications for SZE cluster 
surveys. 


2. THEORETICAL FRAMEWORK 
2 . 1 . Thermal Sunyaev-Zel’dovich Effect 

The thermal SZE is a distortion in the CMB spectrum 
produced by the inverse Compton scattering of CMB 
photons off free electrons in dense structures such as clus¬ 
ters of galaxies. For a given line of sight, a change in the 
CMB specific intensity caused by the thermal SZE at a 
frequency v is given by A/„// C mb = fv{x)g„(x)y. The 
dimensionless comptonization parameter y is defined as, 

y = kB(JT f n e {l)T e (l)dl. (1) 

m e c z J 

where n e and T e are the number density and temper¬ 
ature of electrons, m e is the electron rest mass, c is 
the speed of light, and a B is the Thomson cross-section. 
There are several frequency dependent factors, including 
dcMB = 2 hv 3 /c 2 {e x - l) -1 , /„( x) = [x(e x + l)/(e x - 1) - 
4](1 + <5sze {x,T e )) and g v (x) = x 4 e x /(e x - l) 2 , where 


Jsze (x,T e ) is the frequency dependent relativistic cor¬ 
rection and x = hv/ksTcMB- The corresponding change 
in the CMB temperature is given by AT v /Tcmb = 
f„(x)y. In the Rayleigh-Jeans limit (n <C 200GHz), 
AX^/Tcmb = -2 y and AI„ = (2k B v 2 /c 2 )AT„. 

Let us now consider the SZE signal arising from a 
cluster located at redshift z. The SZE flux integrated 
over a solid angle of observation dfl is given by A S„ = 
f AI„dfl = I cmb fv(x)g u (x)Y , where Y is the integrated 
Compton-y parameter defined as 



1 ( k B a T \ 

d 2 A (z) \m e c 2 J 


UrTedV, 


( 2 ) 


and dfl = dA/d 2 A {z) is the solid angle of the cluster sub¬ 
tended on the sky, d A (z) is the angular diameter distance 
to the cluster, cL4 is the area of the cluster on the sky, 
and dV is the cluster volume. Because Y depends on 
a distance, we will work with the intrinsic thermal SZE 
signal, defined as 


Y int =Yd 2 A {z)(xf gas MT m . (3) 


Note that the integrated SZE flux is linearly sensitive to 
gas mass M gas = f ga sAI and mass-weighted temperature 
T m , where f gas is the gas mass fraction and M is the 
total cluster mass. 


2.2. Self-Similar scaling relations 

In the absence of cooling and heating pr ocesses, clus - 
ters are expected to scale self-similarly llKaiseil R98 6?). 
The self-similar model predicts that the temperature of 
the gas scales with the cluster mass as 

M (xT 3 / 2 E-\z) (4) 

where M = 47rr^A c p cr jt/3 is a halo mass enclosed 
within ta , defined as a radius of spherical volume within 
which the mean density is A r times the critical density 1 
Pcrit) at that redshift llBrvan fc Normanl 11998 1. E(z) 
is the redshift-dependent Hubble parameter, defined as 
H{z) = 100 hE{z) km s -1 Mpc -1 , and it is given by 
E 2 (z) = Om( 1 + z) 3 + Ha for a flat cosmology. 

Inserting Eq.Q]into Eq.[3 we obtain SZE scaling rela¬ 
tions predicted by the self-similar model, 


y int oc 


f gas M 5 / 3 E 2 / 3 {z) 
f eas T 5 / 2 E~\z ) 


( 5 ) 


3. SIMULATIONS 

In this study, we analyze high-resolution cosmological 
simulations of eleven cluster-size systems in the “con¬ 
cordance” flat ACDM model with fl m = 1 — fl A = 0.3, 
fib = 0.04286, h = 0.7 and as = 0.9, where the Hub¬ 
ble constant is defined as 100ft. km s -1 Mpc -1 , and as 
is the power spectrum normalization on 8ft -1 Mpc scale. 
The simulations were done with the Adap tive Re fin ement 
Tree (ART) Wbody+ gasdynamics code llKravtsovll999t 
[Kravtsov et al.lf20028h an Eulerian code that uses adap¬ 
tive refinement in sp ace and time , and (non-adaptive) 
refinement in mass dKlvpin et all 1 200 1) to reach the 
high dynamic range required to resolve galaxy-size ha¬ 
los formed in self-consistent cosmological cluster simula¬ 
tions. 

To set up initial conditions we first ran a low resolu¬ 
tion simulation of 80ft _1 Mpc and 120ft -1 Mpc boxes 
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and selected eleven clusters with mass ranging from 
M 500 C ~ 3.5 x 10 13 to 9 x 10 14 h _1 Mq. Table U lists 
properties of clusters at the present epoch. The pertur¬ 
bation modes in the Lagrangian region corresponding to 
a sphere of several virial radii around each cluster at 2 =0 
was then re-sampled at the initial redshift of Zi =25 for 
three most massive clusters (CL1-3) and Zi = 49 for re¬ 
maining eight clusters in the sample. For the three most 
massive clusters we have resampled radius of 3i? v i r (z=0), 
while for the rest of the clusters the resampling sphere 
had radius of 5i? v ir(z=0), where R w -„{z=Q) is the virial 
radius enclosing overdensities of A v j r = 334 with respect 
to the mean density of the universe at 2=0. During the 
resampling we retained previous large-scale waves intact 
bu t included ad d itiona l small-scale waves, as described 
bv iKlvoin et al.i |200lh . The resampled Lagrangian re¬ 
gion of each cluster was then re-simulated with high dy¬ 
namic range. 

High-resolution simulations were run using 128 3 uni¬ 
form grid and 8 levels of mesh refinement in computa¬ 
tional boxes of 120h -1 Mpc for CL1-CL3 and 80h -1 Mpc 
for CL4-CL11, which corresponds to the dynamic range 
of 128 x 2 8 = 32768 and the peak formal resolution of 
80/32,768 ~ 2.44/i _1 kpc, corresponding to the actual 
resolution of « 2 x 2.44 « 5 h~ x kpc. Only the region 
of ~ 3 — 10/i _1 Mpc around a cluster was adaptively 
refined, the rest of the volume was followed on the uni¬ 
form 128 3 grid. The mass resolution corresponds to the 
effective 512 3 particles in the entire box with a parti¬ 
cle mass to p = 1.07 x lO 9 Af 0 and 3.16 x 10 s Af 0 , or the 
Nyquist wavelength of Any = 0.469h -1 and 0.312/i _1 
comoving megaparsec for CL1-3 and CL4-11, respec¬ 
tively, or 0.018h _1 and O.OO 6 / 1" 1 Mpc in physical units 
at initial redshift of the simulations. The dark mat¬ 
ter particle mass in the region around the cluster is 
therefore (1 — /b)m p = 9.1 x 10 8 h -1 M 0 for CL1-3 and 
2.7 x 10 8 /i _1 M 0 for CL4-CL11 (where fb = = 

0.1429), while other regions were simulated with lower 
mass resolution. 

As the zeroth-level fixed grid consisted of only 128 3 
cells, we started the simulation already pre-refined to the 
2nd level (l = 0,1,2) in the high-resolution Lagrangian 
regions of clusters. This is done to ensure that the cell 
size is equal to the mean interparticle separation and all 
fluctuations present in the initial conditions are evolved 
properly. During the simulation, refinements were al¬ 
lowed to the maximum l = 8 level and refinement cri¬ 
teria were based on the local mass of DM and gas in 
each cell. The logic is to keep the mass per cell approx¬ 
imately constant so that the refinements are introduced 
to follow collapse of matter in a quasi-Lagrangian fash¬ 
ion. For the DM, we refine the cell if it contains more 
than two dark matter particles of the highest mass reso¬ 
lution specie. For gas, we allow the mesh refinement, if 
the cell contains gas mass larger than four times DM par¬ 
ticle mass scaled by the baryon fraction. In other words, 
the mesh is refined if the cell contains DM mass larger 
than 2(1 — fb)m p or gas mass larger than = Af b m p . We 
analyze clusters at the present-day epoch as well as their 
progenitors at higher redshifts. 

We repeated each cluster simulation with and with¬ 
out radiative cooling and star formation. The first set 
of “adiabatic” simulations have included only the stan¬ 


dard gasdynamics for the baryonic component without 
radiative cooling and star formation. The second set 
of simulations included gasdynamics and several phys¬ 
ical processes critical to various aspects of galaxy forma¬ 
tion: radiative cooling, star formation, metal enrichment 
and thermal feedback due to the supernovae type II and 
type la, self-consistent advection of metals, metallicity 
dependent radiative cooling an d UV heating due to cos - 
mological ionizing background llHaardt fc Madatfll996 ). 
We will use labels ’AD’ and ’CSF’ for the adiabatic sim¬ 
ulations and the simulations with gas cooling and star 
formation, respectively. 

Cooling and heating rates take into account Comp¬ 
ton heating and cooling of plasma, UV heating, atomic 
and molecular cooling and are tabulated for the temper¬ 
ature range 10 2 < T < 10° K and a grid of metallicities, 
and UV intensities using the Cloudy code (ver. 96b4, 
iFerland et alJlTflOS ). The Cloudy cooling and heating 
rates take into account metallicity of the gas, which is 
calculated self-consistently in the simulation, so that the 
local cooling rates depend on the local metallicity of the 
gas. Star formation in these simulations was d one us- 
ingan observationally-motivated recipe fe.g.. iKennicuttl 
11998): = p i a 5 /<*, with f, = 4x 10 9 yrs. The detailed 

description and implementation of the star formation are 
provided in Kravtsov et al. (2005). 

The code also accounts for the stellar feedback on the 
surrounding gas, including injection of energy and heavy 
elements (metals) via stellar winds and supernovae and 
secular mass loss. Once formed, each stellar particle is 
treated as a single-age stellar population and its feed¬ 
back on the surrounding gas is implemented accordingly. 
More specifically, in the simulations analyzed here, we 
assumed that a stella_r_ i nitial mass fu nction (IMF) is 
described by the iMiller fc Scald (I1979H functional form 
with stellar masses in the range 0.1 — 100 Mg. All stars 
more massive than AT* > 8 Mg deposit 2 x 10 51 ergs of 
thermal energy in their parent cell 1 and a fraction fz = 
min(0.2, 0.01 AT* — 0.06) of their ma ss as metals, w hich 
crudely approximates the results of lWooslev fc Weaved 
{HU). In addition, stellar particles return a fraction of 
their mass and metals to the surrounding gas at a secular 
rate mi oss = m* Cn(t —thirt.h+Tn )~ 4 with Cq — 0.05 and 
To = 5 Myr 1.Tungwiert et alJ12001 ). The code also ac¬ 
counts for the SNIa feedback assuming a rate that slowly 
increases with time and broadly peaks at the population 
age of 1 Gyr. We assume that a fraction of 1.5 x 10 -2 
of mass in stars between 3 and 8 M 0 explodes as SNIa 
over the entire population history and each SNIa dumps 
2 x 10 51 ergs of thermal energy and ejects 1.3 M 0 of met¬ 
als into the parent cell. For the assumed IMF, 75 SNII 
(instantly) and 11 SNIa (over several billion years) are 
produced by a 10 4 M 0 stellar particle. 

High-resolution 2 and inclusion of various physical pro¬ 
cesses are critical for assessing the effects of galaxy forma¬ 
tion on_ 2 ro 2 erLesand_evolution of the ICM. For exam¬ 
ple, iKravtsov et al.l li2005l) showed that gas cooling and 
star formation significantly suppress gas mass fractions 

1 No delay of cooling was introduced in these cells after SN 
energy release. 

2 The mass and spatial resolution are high enough to resolve 
galact ic-size halos whose mass i s as small as 10“ 4 of the host cluster 
mass tNagai fc Kravtsovl l2003l . 
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TABLE 1 

Simulated cluster sample of the CSF run at z=0 


Name 

Rsooc 

/i _1 Mpc 

M gas 

lvl 500c 

h- 1 10 13 M o 

M tot 

m 500c 

r 1 io 14 M 0 

(Xspec) 

keV 

CLl 

1.160 

8.19 

9.08 

10.0 

CL2 

0.976 

5.17 

5.39 

8.1 

CL3 

0.711 

1.92 

2.09 

4.1 

CL4 

0.609 

1.06 

1.31 

2.7 

CL5 

0.661 

1.38 

1.68 

3.9 

CL6 

0.624 

1.22 

1.41 

3.2 

CL7 

0.522 

0.74 

0.82 

1.7 

CL8 

0.487 

0.43 

0.67 

2.2 

CL9 

0.537 

0.78 

0.90 

2.2 

CLIO 

0.509 

0.62 

0.77 

2.1 

CL11 

0.391 

0.27 

0.35 

1.2 


in clusters. Interestingly, the results of numerical simu¬ 
lations agree qu ite well with recent X -ray and SZE mea¬ 
surements (e.g., LaRoa ue et al.ll200fl Kravtsov et al., in 
preparation). Note also that the level of entropy and 
metal abundance in these simulations also compare well 
with recent X-ray observations. Because the SZE signal 
depends linearly on the gas mass fraction, these results 
indicate that the impact of galaxy formation on the SZE 
flux may also be significant. 

Throughout this paper we use estimates of Y int , mass 
and other cluster observables within different commonly 
used radii, defined by the total matter overdensity they 
enclose. We will use radii r 25 och r soo, ^200 enclosing over¬ 
densities of A c = 2500, 500, and 200 with respect to the 
critical density, p cr it, as well as radii riso and r v ; r enclos¬ 
ing overdensities of A m = 180 and A v ; r with respect to 
the mean density of the universe. The latter is equal to 
A v i r ~ 334 at z=0 and « 200 at z=l for the cosmol¬ 
ogy adopted in our simulations. The virial radius and 
masses of clusters within r$ooc for the CSF run at z =0 
are given in Tabled For reference, we also give a spec¬ 
tral X-ray temperature, (T spec ), of individual clusters ex¬ 
tracted from mock Chandra analysis (Nagai, Vikhlinin & 
Kravtsov 2006, in preparation), and it is equivalent to a 
single-temperature fit to the cluster spectrum extracted 
from the radial range 0.15 < r/r 50 oc < 10. 


4. RESULTS 

In this section, we study the impact of galaxy forma¬ 
tion on the SZE scaling relations using cosmological sim¬ 
ulations of cluster formation. Our general strategy is to 
assess their effects on cluster observables and scaling re¬ 
lations by comparing two sets of simulations with and 
without processes important for galaxy formation. 

4.1. Pressure Profiles 

Since the thermal SZE is sensitive to the pressure of the 
ICM, we will start with the analysis of pressure profiles of 
clusters. Figure ^ shows a comparison of the mean pres¬ 
sure profiles of clusters simulated in the adiabatic regime 
and those with gas cooling and star formation at z=0 and 


z= 1. To compute the mean profiles, we first normalize 
differential and cumulative pressure profiles of individ¬ 
ual clusters using the mean integrated pressure within a 
sphere of r $ooo denoted as -P 5000 and plot them as a func¬ 
tion of the cluster-centric radius in units of r^ooc- The 
mean profiles are computed by averaging the normalized 
profiles over the entire cluster sample. To highlight the 
impact of gas cooling and star formation on the pressure 
profiles, we normalize the pressure profile of the CSF run 
using P 500 C of the AD run for each cluster. The shaded 
bands show 1 <t rms scatter around the mean profile of 
the CSF runs. The mean and scatter of the profiles are 
computed for a logarithm of pressure. 

Fig. □ illustrates that clusters ranging in an order of 
magnitude in mass exhibit remarkably similar pressure 
profiles at all radii. The differential pressure profile de¬ 
clines by nearly four orders of magnitude from the cen¬ 
ter to rsooc, and the scatter is small (~ 20%) through¬ 
out clusters. In the outskirts, the radial profile of gas 
pressure generally falls faster than that of gas or dark 
matter density. This is because pressure is a product of 
gas density and gas mass-weighted temperature for the 
ideal gas, and the latter also declines with radius (see 
iMotl et all 1 20051 iVikhlinin et al. I [20051) . For example, a 
typical power-law slope of the differential pressure pro¬ 
file is ~ —3.5 at r^ooc, which is steeper than a typical 
slope of the density profile of » — 3 see n in cosmological 
simulations llNavarro et alJlToOfl 1199 71) as well as recent 
X-ray observations llVilddinin^'t a 1. 1120061) . 

One of the main results from this analysis is that gas 
pressure is suppressed in a region outside r > 0 . 1 ?' 5 oo c 
in the CSF runs compared to the AD runs. The cumu¬ 
lative pressure at r^ooc, for example, is suppressed by 
~ 25% on average, and this effect is relatively constant 
(~ 25 — 35%) in a region outside > 0.2r 5O o c at z=0. 
In the center, gas pressure is significantly more concen¬ 
trated, becau se of strong gas coo ling and contraction of 
dark matter llOnedin et al l l2004l l. Note, however, that 
this high pressured cluster core contributes very little to 
the cumulative pressure at large radii. For example, the 
inner regions within < O.lrsooc contribute less than 2% of 
the total SZE signal within rsooc- Thus, the SZE signal 
integrated out to a sufficiently large radius is insensitive 
to properties of cluster cores. 

Although general trends are similar, the suppression 
of gas pressure varies more strongly with radius at z=1 , 
because a dense, cool core dominates a larger fraction 
of rsooc at high redshift. Nevertheless, we find that the 
suppression of the cumulative pressure at r$ 00 c is com¬ 
parable at z =0 and 1. However, the difference becomes 
more apparent in the inner region; for example, the effect 
at r 25 ooc changes from 35% at z =0 to 20% at z= 1. If a 
similar trend exists in real clusters this implies that red- 
shift evolution of the SZE scaling relations is expected 
to deviate from the prediction of the self-similar model, 
if the SZE flux is computed within the aperture consid¬ 
erably smaller than r$ ooc- 

In order to better understand the suppression of gas 
pressure, we examine the impact of gas cooling and star 
formation on density, temperature and entropy profiles 
of the ICM. Figure [3 shows that the suppression of pres¬ 
sure is accompanied by strong suppression in gas density 
and moderate increase in gas mass-weighted tempera¬ 
ture. An entropy profile provides further insights into 
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Fig. 1 . — Cumulative (top panels) and differential (bottom panels ) pressure profiles for the eleven clusters used in our analysis at 2=0 
(left column) and 2=1 (right column). The dashed and solid lines show mean profiles in the AD and CSF simulations averaged over all 
clusters, respectively. Note that all profiles are normalized to the mean cumulative pressure at rsooc in the AD runs, Psooc, at each epoch. 
The bottom panel of each figure shows the fractional deviation of the mean profile of the CSF run relative to the mean profile of the AD 
run. In all panels, the shaded bands indicate 1 o rms scatter around the mean. The vertical arrows in the top panels indicate the radii 
enclosing overdensities of 2500 and 200 with respect to the critical density at each epoch. 


these phenomena, because it uniquely characterizes ther¬ 
modynamic properties of the ICM, given a shape of its 
confining potential well (see Voit et al. 2002 and refer¬ 
ences therein). Fig.[2shows that the entropy of the ICM, 

defined as S = T/nl^ 3 , is enhanced from 0.05 x rsooc out 
to 2 x r 5 ooc- The effect peaks around 0.3r5oo c , and the 
entropy is enhanced by about 40% at rsooc- In our sim¬ 
ulations, gas cooling and star formation can increase the 
entropy of the gas in one of two ways: (1) the level of 
entropy is elevated throughout as gas co oling and star 
formation proceed in the cluster center llVoit fc Brvanl 
l200l 3 or (2) direct heating of the ICM by energy in- 

3 Gas cooling and star formation remove low entropy gas and 


jection from supernova explosions. The former is likely 
a dominant mechanism that shapes the overall entropy 
distribution, while the second process is important for 
regulating a rate of gas cooling in high density regions 
(e.g., a vicinity of a central galaxy and cluster galaxies), 
where stars are forming with sufficiently high rates. 

To examine a mass dependence of the effects, we split 
the sample in half and show the mean profiles of six most 
massive (Tx > 2keV) and five least massive (Tx < 2keV) 

cause an inflow of high entropy gas from outside to maintain the 
stability of the cluster core. In this way, the low-entropy gas in 
the inner region is gradually replenished with the high-entropy gas 
flowing in from outside, and the level of entropy is gradually el¬ 
evated throughout the cluster as gas cooling and star formation 
proceed in the cluster core. 
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Fig. 2. — The impact of galaxy formation on pressure, gas den¬ 
sity, temperature and entropy profiles (from top to bottom) of the 
simulated clusters at z=0. The solid lines show fractional changes 
in the mean differential radial profiles in the CSF runs averaged 
over eleven clusters relative to the mean profiles in the AD runs, 
indicated by lines at zero. The dotted and dashed lines show the 
same for subsamples of six most and five least massive clusters, 
respectively. Note that the mean profile of each quantity is com¬ 
puted by first normalizing a profile of each cluster by the value of 
the cumulative profile at rsooc and averaging the normalized profile 
over all clusters. 


clusters in Fig. [3 Note that the former includes clusters 
dominated by bremsstrahlung or those that are about 
equal in the importance of bremsstrahlung and line cool¬ 
ing, while the latter includes clusters dominated by line 
emission. Fig. |3 shows that a shape and normalization 
of the pressure profiles exhibit the smallest systematic 
trend with mass. In the outskirts (r > 0.3r5oo c ) of 
the clusters, there is a similarly small systematic trend 
with mass for gas density, temperature and entropy pro¬ 
files. In the inner regions, on the other hand, we find 
a slightly more pronounced mass dependence in all pro¬ 
files. Among them, the pressure profiles show the small¬ 
est trend with mass. In less massive clusters, changes in 
gas density and temperature become larger but with op¬ 


posite signs; therefore, these effects cancel each other to 
give the smallest mass dependence on pressure. In con¬ 
trast, the effects add constructively to give the largest 
mass dependence on entropy. 

The main conclusions from the profile analysis are that 
gas cooling and star formation suppress the amplitude of 
gas pressure and therefore the SZE signal, and the mag¬ 
nitude of the effects depends rather weakly on the clus¬ 
ter mass. These results suggest that processes of galaxy 
formation modify normalization of the SZE flux-mass re¬ 
lation, but its impact on slope is expected to be small. 
This is a main topic of the discussion in the next section. 

4.2. The SZE scaling relations 

The SZE flux-total mass (SZE-M) relation is a relation 
most directly relevant for cosmological application. Fig¬ 
ure 0 shows a comparison of the SZE-total mass scaling 
relation at r 5 oo c and its redshift evolution for a sample of 
eleven clusters in the AD and CSF runs. At each epoch, 
we performed fits to the sample of simulated clusters us¬ 
ing the simple power-law relation, 

Y = A "* 10 ~ 6 (lo^FT tjj) . (6) 

where Au is the normalization at 10 14 H~ 1 Mq in units of 
10 -6 and is the slope. In practice, we fit a straight 
line to the log(Y)-log(M) relation by minimizing y 2 . Ta¬ 
ble 12141 lists the best-fit normalization and slope mea¬ 
sured at different radii, including nsom, r v i r , r 2 ooc, rsooc 
and r 25 ooc- The values in the adiabatic simulations are 
marked ’ad’, while those in the simulations with gas cool¬ 
ing and star formation are marked ’csf’ (e.g., The 

best-fit parameters in these tables and the right panel of 
Fig.0 are obtained by fitting for both normalization and 
slope simultaneously. To highlight a deviation from the 
self-similar slope, the best-fit relations shown in the left 
panel of Fig. 0 are obtained by fixing the slope to the 
predicated value of 5/3. Note that fitting for the slope 
changes the best-fit normalization by no more than 5% 
at both 2=0 and 2 = 1 . 

In the right panel of Fig. [3J errorbars indicate 2tr con¬ 
fidence region of the best-fit slope and normalization in 
the CSF simulations at each epoch. Since the redshift 
evolution of normalization is measured relative to the re¬ 
lation at 2 = 0 , errorbars at 2 >0 include uncertainties in 
the best-fit normalizations at each epoch and 2=0 added 
in quadrature. Note also that the size of the errorbars 
are comparable for the AD simulations. Errorbars shown 
in Fig. M and Fig. Q are also computed in the same way. 

The left panel shows that the SZE signal integrated 
within a sphere of r$ 00 c correlates very strongly with the 
enclosed cluster mass at both 2=0 and 1. The rms scatter 
is about 10-15% in both AD and CSF runs. Moreover, we 
find that the slope of the AD run is in a very good agree¬ 
ment with the predicted slope of the self-similar model. 
The best-fit slope of the CSF run is also not very dif¬ 
ferent from the predicted slope of 5/3, but there is an 
indication that the slope may be systematically steeper 
by 1$ 0.1. Unfortunately, however, the current sample 
size is too small to rule out the prediction of the self¬ 
similar model with any statistical significance. A larger 
sample of simulated clusters is clearly needed to assess 
this effect in future studies. 
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Fig. 3.- The SZE flux-total mass relations for eleven clusters in the AD (solid symbols) and CSF (open symbols) simulations. Left-top 
panel: the relation between the integrated compton-Y vs. total mass enclosed within a sphere of radius rsooc at z=0 (bottom lines and 
points) and z =1 (top lines and points). For clarity, the z=l relations are shifted upward by a factor of ten. The solid circles and solid 
triangles indicate the AD runs at z =0 and 1, and the open squares and open stars are the CSF runs at z = 0 and 1, respectively. At 
each epoch, the dotted and dashed lines show the best-fit relations for the AD and CSF runs with a slope fixed at the self-similar slope 
of 5/3. Left-bottom panel: the fractional deviation of individual clusters from the best-fit relation of the AD runs at each epoch. The 
dashed line indicates the best-fit relation of the CSF runs at z= 0. Right panel: the redshift evolution of slope (top) and normalization 
(middle) of the SZE flux-total mass relation between z=2 and the present day. The solid circles and open squares indicate the AD and 
CSF runs, respectively. The dashed lines indicate the slope of 5/3 (top) and the evolution in normalization predicted by the self-similar 
model (middle). The errorbars indicate 2 a confidence region of the best-fit slope and normalization in the CSF runs at each epoch. The 
bottom panel shows a fractional deviation of the best-fit normalization from the self-similar prediction as a function of time. 


While the tightness and slope of the relations are rel¬ 
atively unaffected, gas cooling and star formation sig¬ 
nificantly modify normalization of the SZE-total mass 
relation. Compared to the AD run, the normalization of 
the CSF run is lower by 41, 34 and 27% at r 2500 c, ^ 500 c 
and r 180m . 

The right panel shows that the redshift evolution of 
the SZE-total mass relation is consistent with the pre¬ 
dictions of the self-similar model in the entire redshift 
range considered here. First, the slope is constant with 
redshift. Second, the redshift evolution of the normaliza¬ 
tion is in a good agreement with the predicted evolution. 
We also find similar results at other radii considered in 
this paper. 

To better understand the impact of galaxy formation 
on the SZE observable-mass relation, we also examine 
their effects on f gas and T m in Figure^ Compared to the 
AD run, f gas is lower by 38% in the CSF run; however, its 
effect on the SZE flux is offset partially by an increase in 
T m by 9%. The net effect is a reduction of the SZE signal 
by 32%, which falls slightly short of a 34% change in the 
normalization of the SZE-M relation. This illustrates 
that changes in f gas and T m alone do not fully account 
for the change in normalization. 

In addition, gas cooling and star formation slightly 
modify the total cluster mass. Figure 0 illustrates that 
inclusion of these processes cause an increase in M tot at 
r 500 c by about 6% on average. The effect becomes larger 
at a smaller radii, and it is about 10% at r 25 ooc- This 


effect causes shifts in the SZE-M relation in both x- and 
y-directions by the same amount. However, since the 
slope of the relation is ~5/3, the 6% shift in the positive 
x-axis translates into a 8% shift in the negative y-axis. 
The net effect is therefore a negative 3% change in nor¬ 
malization due to the change in M tot . Putting them all 
together, we find that the 34% decrease in normalization 
of the SZE-M relation at r 50 oc is due to combined effects 
of the 38% decrease in f gas , the 9% increase in T m and 
the net 3% decrease in normalization due to the change 
in M t ot- Note also that inclusion of gas cooling and star 
formation lowers the integrated SZE signal ( Y mt ) and 
the gas mass (M g ) within r 50 oc by 29% and 35%, respec¬ 
tively. 

Finally, we check the redshift evolution of f gas and 
T m . The self-similar model predicts that f gas is constant 
with time, while T m of a given clutter mass evolves with 
time according to Eq. 0] To examine deviations from the 
predicted evolution, Figure El plots the deviations of the 
best-fit normalization of the f gas — M and the T m — M 
relations at M^ooc = 10 14 /i _1 M 0 from the self-similar 
prediction as a function of redshift for the AD and CSF 
runs. Fig. mi shows that f gas of a given mass in the AD 
runs is consistent with the predicted evolution of the self¬ 
similar model, while gas cooling and star formation cause 
deviations, which increase toward higher redshifts. In 
the CSF runs, f gas is higher by about 12% at z=l and 
17% at z =2 compared to their values at z= 0. T m , on the 
other hand, is lower by an almost equal amount at higher 
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Fig. 4. — Fractional changes in gas mass fraction (top) and mass- 
weighted temperature ( bottom ) enclosed within a sphere of radius 
7*500c between the AD and CSF simulations as a function of cluster 
mass. The solid circles and solid triangles are the AD runs at 
z =0 and 1, and the open squares and open stars are the CSF runs 
at these epochs, respectively. The dashed line indicates a mean 
fractional change in f ga s or T m in the CSF run at z=0 compared 
to that of the AD run, indicated by the dotted line at zero. 
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Fig. 5. — Effects of gas cooling and star formation on total mass 
profiles of simulated clusters at z=0. The solid line shows a change 
in the total mass profiles in the CSF simulations compared to that 
in the AD simulations averaged over ten clusters. For each cluster, 
we compute a fractional change in the total mass profiles of the 
CSF and AD runs, (Mcsf — Ma)/Ma, at the same physical radius, 
where Mcsf and Ma are the total mass profiles of the CSF and 
AD runs. We then plot the fractional change between two runs as 
a function of radius in units of r*500c and average over ten clusters, 
while excluding one cluster experiencing a major merger at z=0. 


redsliifts in both the AD and CSF runs. Since the SZE 
signal is linearly proportional to both of these quantities, 
the evolution in f gas is canceled almost exactly by the 
evolution in T m . Note also that we did not find any sys¬ 
tematic evolution in the impact of gas cooling and star 
formation on the total mass with redshift. This explains 
why the SZE signal of a given mass shows very little evo- 


Fig. 6 .— The evolution of gas mass fraction (top) and mass- 
weighted temperature ( bottom ) at M500C = 1 0 14 h ~ 1 A 7 q . The frac¬ 
tional deviation of the best-fit normalization from the self-similar 
prediction, indicated by the dashed line. The errorbars indicate 2cr 
statistical uncertainties in the CSF runs at each epoch. 


lution beyond the evolution predicted by the self-similar 
model. 

So far, we have focused on the SZE-total mass relation 
because of its direct relevance for cosmological studies. 
Unfortunately, however, the cluster mass is not directly 
observable, and making unbiased measurements of the 
SZE-M relation is difficult in practice. Therefore, it is 
also useful to study relations between the SZE signal and 
cluster properties that can be measured more reliably 
from observations, as they can provide useful direct tests 
of simulations. 

Here, we consider two such quantities: gas mass (M g ) 
and mass-weighted temperature ( T m ) of clusters. Fig¬ 
ure [3 shows SZE -M g and SZE-T m relations at r^ooc and 
evolution of their slope and normalization with redshift. 
At each timestep, we perform a fit to the entire sample of 
simulated clusters using the following simple power-law 
relations, 


Y = A 13 x 1CT 6 

Y = A 5 x 1CT 5 I 


( M g 
v 10 13 h^M Q 
T \“ T 

- L m \ 

5keV ) 


a Mg 


(7) 

( 8 ) 


where A 13 is the normalization at 10 13 h 1 Mq in units of 
10 -6 and aM g is the slope of the SZE-M g relation. Sim¬ 
ilarly, A 5 is the normalization at 5 keV in units of 10 -5 
and ap is the slope of the SZE-T m relation. Table EID list 
the best-fit normalizations and slopes measured at differ¬ 
ent radii for a z =0 sample. The errorbars indicate a 2cr 
confidence region of the best-fit slope and normalization 
in the CSF runs at each epoch. 

Fig. [3 shows that both the SZE -M g and the SZE-T m 
relations are as equally tight as the SZE-M relation. The 
best-fit slopes of the AD runs are in a very good agree- 
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Fig. 7.— Top panels: the SZE flux-gas mass relations. The left panel shows the scaling relations at z =0 (bottom lines and points ) and 1 
(top lines and points). The right panel shows the redshift evolution in slope and normalization between z =2 and 0. The errorbars indicate 
2rr confidence region of the best-fit slope and normalization in the CSF runs at each epoch. Point and line types are the same as in Fig. FTI 
Bottom panels: the relation between the SZE flux and the mass-weighted temperature. 


ment with the self-similar slope. The best-fit slopes of 
the CSF runs, on the other hand, are marginally consis¬ 
tent with the predicted slope, and there are indications 
that the slopes in the CSF run may be slightly smaller 
(larger) for the SZE- M g (SZE- T m ) relations. We also find 
that the slope is constant with time in both the AD and 
CSF runs. 

Existence of tight relations between the SZE observ¬ 
able and cluster masses and temperature is encouraging 
news for future cosmological studies. In fact, the tight¬ 
ness of these relations is rather remarkable, because the 
samples analyzed here include clusters in a wide range 
of dynamical states. This indicates that the integrated 
thermal SZE signal is not very sensitive to dynamical 
states of clusters, making the thermal SZE an excellent 
proxy of cluster mass. The reason that the thermal SZE 
works so well is as follows. Fundamentally, the ICM is 
trapped in the external potential of the dark matter and 


the gas must remain in approximate hydrostatic equi¬ 
librium in the dark matter potential provided that the 
perturbations are small compared to the energy scale of 
the cluster itself. The pressure profile is therefore deter¬ 
mined by approximate force balance against gravity, and 
the thermal SZE measures the integral of that profile. 

While the tightness and slopes of these relations are rel¬ 
atively unaffected, gas cooling and star formation have a 
large impact on their normalization. Fig. 0 shows that 
the normalization of the CSF run is higher by 49% in the 
SZE-M g relation at 7'5oo c , while it is lower by 39% in the 
SZE-T m relation. Note that the change in the normal¬ 
ization of the SZE -M g relation has an opposite sign from 
the change in the SZE-M relation, because gas cooling 
and star formation significantly reduce cluster gas mass. 
Recall that the inclusion of these physical processes re¬ 
duces the SZE signal by 29% and M g by 35% at r 50 o c ■ 
For the slope of 5/3, a decrease in M g by 35% in the x- 
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TABLE 2 

Normalization and slope of the SZE-Total Mass 

RELATIONS AT Z =0 


TABLE 4 

Normalization and slope of the SZE-Temperature 
RELATIONS AT z= 0 


A 

AE id 

14 

ad 

a M 

4 cs f 

^14 

rv csf 

180m 

2 ^Q+ 0 - 31 
0.28 

1.63±0.09 

1 CO - ! - ®-! 8 
i.QO_0 17 

1.68±0.07 

vir 

oii +0.38 

0.ii_o.34 

1.61±0.09 

9 orj+0.22 

z.zU—0.20 

1.67±0.08 

200c 

q 'TO+0.42 
' °-0.38 

1.62±0.10 

2 5R+ 0 - 24 

z.ou—0.22 

1.70±0.08 

500c 

4 99"*"5'1!2 

1.66±0.09 

q qq+0.26 

O.Zb»_ 0> 24 

1.73±0.08 

2500c 

q QO + 0.65 

8.82_o 60 

1.68±0.07 

c 99+0.45 
O.z^_o 41 

1.77±0.09 




TABLE 3 



Normalization and slope of 

the SZE-Gas Mass 


RELATIONS AT 

z =0 


A 

4 ad 
^13 

n ad 

a Mg 

4 cs f 

a 13 

“Ms 

180m 

1 cc+0.23 
l.OO — 0 20 

1.63±0.09 

2 1Q+ 0 - 22 
Z, -- Ly — 0.20 

1.60±0.07 

vir 

l-86io.22 

1.61±0.09 

2 61“*" 0,23 
z ' ol -0.21 

1.60±0.07 

200c 

9 9Q+0.22 
-0.20 

1.62±0.07 

3 24”*” 0,21 
-0.20 

1.60±0.06 

500c 

3-08 ±8;1I 

1.65±0.08 

4 59+ 0 28 
<±.oy_0 27 

1.60±0.06 

2500c 

6-3it8:£ 

1.64±0.07 

q 97+0.60 
-0.56 

1.58±0.05 


direction corresponds to an increase in normalization by 
a factor of two. Combined with a reduction of the SZE 
signal, the net effect is an increase in normalization by 
49% at r 50 o c . Note that the effects become larger in the 
inner region of the clusters. Similar accounting works for 
the SZE-T m relation. 

In contrast to the SZE-M relation, gas cooling and star 
formation modify the redshift evolution of normalization 
in the SZE -M g and the SZE-T m relations. In the CSF 
run, the normalization of the SZE -M g relation is system¬ 
atically lower by ~20% at z =1 than the prediction of the 
self-similar evolution model, while the SZE-T m relation 
is higher by w40% at z= 1. The departure from the pre¬ 
dicted evolution is mainly due to the evolution in f gas 
and T m , not the SZE signal, as discussed above (see also 
Fig. [HI • Note also that the redshift evolution of the AD 
simulations is consistent with the self-similar model at 
all redshift. 


4.3. Comparisons to Previous Work 

We now compare the results of our sim ul ations 
to previous studies in literature ffWhite et al 3 1200 ^ 
Ida Silva et all120041: lMot.1 et ai~ll2005|b So far, all of the 
previous studies have focused on simulating a large num¬ 
ber of clusters. While these simulations provide good 
statistics, the resolution is inevitably limited to capture 
relevant physical processes. The present work is on the 


A 

A|d 

C+ 



180m 

18.24 jig'®® 

2.52±0.18 

12.06 t 3 2 f 7 

2.60±0.21 

vir 

14 11 +2-80 
1^-11_2.34 

2.52±0.15 


2.63±0.19 

200c 


2.50±0.13 

6.37i l:ll 

2.62±0.17 

500c 

5 Q7+0-84 
0>y '-0.74 

2.45±0.12 

q fjq+0.62 
O.UO-0 53 

2.61±0.17 

2500c 

2 1Q+ 0 - 37 
z,±u -0.32 

2.50±0.17 

1 1 9+0-18 
i-iy_o.i6 

2.73±0.19 


other end of the spectrum: the resolution is high, but the 
statistic is low. In this sense, our work is complimentary 
to the previous studies. 

Using SPH simulations, Ida Silva et ali ( I2004fl studied 
the SZE scaling relations in adiabatic simulations and 
re-simulations in which gas was allowed to cool or was 
pre-heated. These studies showed that both strong gas 
cooling and preheating modify the slope of the SZE-M 
relation. They find that the slope in their cooling run 
(qm = 1.79) was steeper than the slope of their adiabatic 
run (a at = 1.69) by about 0.1. More recently. iMotl et all 
( 2005') used a large sample of clusters simulated using an 
Eulerian AMR code and studied the SZE-M relation in 
the adiabatic run and three resimulations in which they 
added gas cooling, star formation and feedback one at a 
time. These studies showed that the slope in the simu¬ 
lations that include both gas cooling and star formation 
is very similar to the slope of the adiabatic run. Similar 
results are obtained using simulations performed with 
the entropy conserving G adget code with cooling, star 
formation and feedback {White et a,1.ll2f)fl3 i. Results of 
our simula tions are al so consistent with these findings. 
Moreover, IMotl et abl ( 2005 ) showed that the slope in 
the cooling only run is steeper than the slope in the adi- 
ab atic run b y abo ut 0.1 , consistent with results reported 
bv lda Silva, et al ](l200l. 

In addition. Ida Silva et al 1 d200l investigated the red¬ 
shift evolution of the slope and scatter. They showed 
that the slope and normalization of the SZE-M relation 
evolve with redshift according to the self-similar model 
out to z = 2 in th eir adiabatic and cooling only runs. 
iMotl et ah J20051) also found similar results in all four 
sets of their simulations up to ^=1.5. We reached the 
same conclusions using our high-resolution simulations. 

In summary, the scatter, slope and redshift evolution 
of the SZE observable-mass relation are generally insen¬ 
sitive to details of cluster gas physics. These results ap¬ 
pear to be very robust. In fact, the agreement among 
different studies is rather remarkable, because these sim¬ 
ulations were carried out using very different numeri¬ 
cal techniques, resolution and implementation of various 
physical processes incorporated in simulations. 

Despite the robustness of the results discussed so far, 
the impact of galaxy formation on normalization of the 
SZE scaling relations is not yet well understood, be¬ 
cause the effect is closely related to the cold baryon frac¬ 
tions, which have not yet converged among different sim¬ 
ulations. Ida Silva et al. (1200 4th for example, find that 
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Fig. 8. — Comparisons of the SZ flux-gas mass relations in the 
simulations and observations. The data points with errorbars show 
the BIMA/OVRO S Z+Chandra X- ray cluster observations for a 
sample of 36 clusters JLaRoauei2Q05l) . The SZ flux is the integrated 
flux within the 2D projected aperture of R 2500 c, while the gas mass 
is the enclosed gas mass within the sphere of r* 2500 c- The SZ flux 
is also corrected for the redshift evolution assuming the self-similar 
evolution model and the ACDM cosmology. The solid circles and 
open triangles show the eleven simulated clusters in the adiabatic 
and cooling and star formation runs, respectively, at z=0. The 
dotted and dashed lines are the best-fit relations at z=0 of each set 
of simulation. The figure shows that the simulations that include 
gas cooling and star formation are in reasonably good agreement 
with the data, while the adiabatic simulations are inconsistent with 
the observed relation. 


the normalization of the SZE-M relation in their cooling 
only run is lower by 20% and 9% at M 2 oo c = 10 14 and 
5 x 10 14 /i” 1 Mq com pared to the adi abatic run, respec¬ 
tively. Moreover, I White et al.l (' 2002 1 find even smaller 
effects in their simulation, because of the efficient en¬ 
ergy feedback model used to regulate conversion of hot 
gas into stars llSpringel fc Hernauistl 12002 1. Using our 
high-resolution simulations, we report that processes of 
galaxy formation lower the normalization of the SZE-M 
relation by 41, 34 and 32% at r 2 500 c> ^ 500 c and r 2 ooc, re¬ 
spectively, almost uniformly throughout the cluster mass 
range. This is one of the largest impact of galaxy forma¬ 
tion on the SZE scaling relations reported to date. 

4.4. Comparisons to Observations 

To gauge how well the current simulations compare to 
data, we compare the scaling relations in our simulations 
to recent observational results. Figure 0 shows a com¬ 
parison of the SZE flux-gas mass relation from our sim¬ 
ulations to recent observational results based on a sam¬ 
ple of 36 clusters between 0=0.14 and 0=0.89, observed 
with the SZE BIMA /OVRO interfer ometers and Chan¬ 
dra X-ray telescope llLaRoauel 1200511) . For comparison, 
we use the SZE flux-gas mass relation, rather than the 
SZE-total mass relation, because gas mass can be mea- 
sure d more accurate ly using Chandra X-ray observations 
(see lLaRoaudll2005 !. The SZE flux is an integrated flux 


within a 2D projected aperture of I? 2500 c while the gas 
mass is the enclosed gas mass within a sphere of r 2 5 ooc- 
These quantities are computed using the best-fit parame¬ 
ters of the isothermal 0 mod el, fi t joint ly to the SZE and 
the Chandra data fsee lLaRoaudl2005l for details). The 
observed SZE flux at a given redshift is also corrected for 
the redshift evolution assuming the canonical self-similar 
model and the ACDM cosmology. 

In the simulations, the SZE flux is computed by pro¬ 
jecting a sphere of 3xrsoo c centered around the minimum 
of a cluster potential. Note that hot gas in the cluster 
outskirt makes non-negligible contribution to the inte¬ 
grated SZE signal, if the projected aperture is a small 
fraction of the cluster virial radius. For a projected 
radius of R 25000 approximately 35% of the SZE signal 
arises from a region outside a sphere of r 2500 c on aver¬ 
age. The projection effect varies by about 20% among 
different clusters, but we did not find systematic varia¬ 
tion of the effect with cluster mass. The contribution of 
the cluster outskirt becomes less significant as we make 
the projected aperture larger. The projection of the hot 
gas associated with dense structures in the foreground 
or background of the cluster is not accounted for in this 
comparison. 

The comparison shows that the adiabatic simulations 
are strongly inconsistent with data. Inclusion of gas cool¬ 
ing and star formation brings simulations in a better 
agreement with observations. However, a discrepancy 
between simulations and observations still remains. More 
specifically, there are indications that the normalization 
is larger, and the slope is shallower for the observed 
relation fsee lLaR,oauel 12005 ). However, the remaining 
discrepancy is at a level of sys tematic uncertainties in 
current measurements (see e.g.. iBenson et alj 1 2004 ). It 
is therefore critical to resolve systematic uncertainties 
among different instruments and to increase the cluster 
sample especially at low-masses to provide a better lever¬ 
age on the slope of the relations. 

5. DISCUSSION & CONCLUSIONS 

We have presented the analysis of the Sunyaev- 
Zel’dovich effect (SZE) scaling relations using high- 
resolution simulations of galaxy clusters formed in a con¬ 
cordance ACDM cosmology. The simulations of eleven 
individual clusters spanning a decade in mass (M 500 C = 
3.5 x 10 13 to 9 x 1O 14 /i _ 1 M0) are performed with the 
shock-capturing Eulerian adaptive mesh refinement N- 
body+gasdynamics ART code. We study the effects of 
gas cooling and star formation on the SZE scaling rela¬ 
tions and their redshift evolution between 0=0 and 0=2 
by comparing two sets of simulations performed with and 
without these processes included. The main results are 
summarized as follows. 

1. The SZE signal integrated to a sufficiently large 
fraction of cluster volume correlates very strongly 
with the enclosed total cluster mass, independent of 
details of gas physics and dynamical state of clus¬ 
ters. The rms scatter of the SZE-total mass relation 
is about 10-15%. 

2. The slope of the relation in the adiabatic run is in 
a very good agreement with the predicted slope of 
the self-similar model in the entire redshift range. 
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We find that the impact of galaxy formation on the 
slope is small (< 0.1). 

3. The redshift evolution of the SZE-total mass re¬ 
lation is consistent with the self-similar model be¬ 
tween z =0 and z= 2: (a) the slope is constant with 
redshift and (b) the normalization evolves with red- 
shift according to the self-similar evolution model. 

4. Gas cooling and star formation significantly modify 
the normalization of the SZE flux-total mass rela¬ 
tion. Inclusion of these physical processes causes a 
decrease in the normalization by 41, 34 and 27% at 
f 2500 cj r^ooc and risom, respectively. The decrease is 
due to a large decrease in gas fraction, which is offset 
somewhat by an increase in mass-weighted temper¬ 
ature. Gas cooling and star formation also cause an 
increase in total cluster mass and hence modify the 
normalization by a few percent. 

5. The integrated SZE signal also correlates strongly 
with gas mass and mass-weighted temperature of 
clusters. The results (l)-(3) apply equally well for 
these relations, except that the redshift evolution of 
the normalization exhibits some deviations from the 
self-similar model, which increases toward higher 
redshifts. Gas cooling and star formation also signif¬ 
icantly modify the normalization of these relations. 

6. The SZE flux-gas mass relation in the simulations 
with gas cooling and star formation is in a bet¬ 
ter good agreement with the observed relation for 
a sample of 3 6 OVRO/BIMA SZE + Chandra X-ray 
observations (iLaRociu IIS than the simulations 
neglecting galaxy formation. 

These results have a number of important implica¬ 
tions for cosmological studies with upcoming SZE clus¬ 
ter surveys. First and foremost, the SZE fluxes of clus¬ 
ters exhibit a remarkable regularity at all redshifts, and 
the SZE signal integrated to a sufficiently large fraction 
of c luster volume is insensitive to merging events (see 
e. g.,lMotl et al]l2005l) or properties of a cluster core (see 
§ S3)- This indicates that the integrated SZE flux is 
an extremely good proxy for cluster mass. Second, the 
slope and the redshift evolution of the SZE scaling rela¬ 
tions are insensitive to details of cluster gas physics, and 
they are well characterized by a simple self-similar clus¬ 
ter model between z =0 and 2. The simplicity of t heir 
redsh i ft evolu t ion implies that the self-calibration (fHul 
12003 ; Maiumdar & Mohil2004f) will be effective. Finally, 
these results appear to be very robust, as the same con¬ 
clusions have been reached using simulations with very 
different numerical techniques, resolution and implemen¬ 
tation of various physical processes incorporated in sim¬ 
ulations (see §m. 

Despite the simplicity of redshift evolution, the nor¬ 
malization of the SZE flux-mass relation is much less un¬ 
derstood, because the effect is closely related to the cold 
baryon fraction, which has not yet converged among dif¬ 
ferent simulations. Using high-resolution cluster simula¬ 
tions, we show that processes of galaxy formation have a 
significant impact on the normalization of the SZE scal¬ 


ing relations. Gas cooling and star formation suppress 
the normalization by ~ 30 — 40%, primarily due to a 
large reduction in cluster gas mass fraction. Interest¬ 
ingly, the SZE scaling relations in these simulations are 
in a reasonably good agreement with recent observations. 
Moreover, the gas mass fractions in these simulations 
compare well with measurements from deep Chandra X- 
ray observations of nearby clusters (Kravtsov et al. in 
preparation). Despite these successes in matching SZE 
and X ray observations, the current cluster simulations 
may still suffer from the ’’overcooling” problem since the 
fraction of baryons in the cold gas and stars within a 
virial radius at z =0 is in a range 0.25-0.35, at least 
a fa ctor of two higher tha n observational measurements 
fsee lKravtsov et al.ll2005l and discussions and references 
therein). 

Given the importance of these issues, further efforts 
from observers and theorists are needed to better under¬ 
stand the SZE scaling relations for cosmological studies. 
Observationally, it is important to increase the sample 
size and the number of low-mass clusters. Comparisons 
of different instruments will also help resolve system- 
atic u ncertain t ies a mong different instruments (see e.g., 
iBenson et al . 3QQl. It is also critical to understand sys¬ 
tematic uncertainties in measurements of cluster mass 
through detailed and extensive comparisons of X-ray, 
SZE and optical observations. With the advent of a 
number of dedicated SZE survey instruments, the ob¬ 
servational situation is expected to improve rapidly. 

Theoretically, it is important to push detailed theo¬ 
retical modelling of cluster formation to investigate roles 
of various physical processes, such as thermal conduction 
and AGN feedback, in shaping the properties of the ICM. 
For cosmological applic ation, it is also cri tical to under¬ 
stand projection effects ifWhite et al.l2002|l and nature o f 
scatter in observable-mass relations ibima, X Hnll2flfl5 ). 
Numerical simulations of cluster formation will likely 
provide important insights into these issues and help as¬ 
sess t he effectiveness o f the self-calibration technique (l]u 
12003 : Maiumdar fc Mohil 12004 1 for future cluster sur¬ 
veys. 
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